DEMO-INP-HEPP-2011-4 
WUB/11-13 
TTK-11-43 

IFJPAN-IV-2011-7 



HELAC-NLO 

G. Bevilacqua a , M. Czakon a , M.V. Garzelli b,c , A. van Hameren d , 
A. Kardos b , C.G. Papadopoulos 6 '*, R. Pittau f , M. Worek s 

a Institut fur Theoretische Teilchenphysik und Kosmologie, RWTH Aachen University, 

D-52056 Aachen, Germany 
b Institute of Physics, University of Debrecen, H-J^OIO Debrecen P.O.Box 105, Hungary 
c University of Nova Gorica, Laboratory for Astroparticle Physics, SI-5000 Nova Gorica, 

Slovenia 

d The H. Niewodniczanski Institute of Nuclear Physics, Polisch Academy of Sciences, 
Radzikowskiego 152, 31-34-2 Cracow, Poland 
^Institute of Nuclear Physics, NCSR Ar\\±6xpiToz, 15310, Athens, Greece 
J Departamento de Fisica Teorica y del Cosmos, Universidad de Granada, E-18011 

Granada, Spain 

9 Fachbereich C, Bergische Universitat Wuppertal, D-^2091 Wuppertal, Germany 



Abstract 

Based on the OPP technique and the HELAC framework, HELAC-1L00P is a 
program capable to numerically evaluate QCD virtual corrections to scatter- 
ing amplitudes. A detailed presentation of the algorithm, along with instruc- 
tions to run the code and benchmark results are given. The program is part 
of the HELAC-NLO framework that allows for a complete evaluation of QCD 
NLO corrections. 
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Nature of problem: 

The evaluation of virtual one-loop amplitudes for multi-particle scattering is a long- 
standing problem (1). In the last years the OPP reduction technique (J2j) opened 
the road for a fully numerical approach based on the evaluation of the one-loop 
amplitude for well-defined values of the loop momentum. 
Solution method: 

By using HELAC (3-5) and CutTools (6), HELAC-1L00P is capable to evaluate QCD 
virtual corrections ((7|). The one- loop n-particle amplitudes are constructed as part 
of the n + 2 tree-order ones, by using the basic recursive algorithm used in HELAC. 
A Les Houches Event (LHE) file is being produced combining the complete in- 
formation from tree- and virtual one-loop contributions. In conjunction with real 
corrections, obtained with the use of HELAC-DIPOLES the full NLO corrections 
can be computed. The program has been successfully used in many applications. 
Additional comments: 



Program obtainable from: http://helac-phegas.web.cern.ch/helac-phegas 
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1. Introduction 



With the advance of LHC experiments, more precise theoretical predic- 
tions will be indispensable. In the last years many groups have been able to 
compute NLO corrections for multi-particle processes a task thought almost 
impossible before 

(0-E3). The NLO revolution {H became possible be- 
cause new reduction techniques have been proposed and implemented, among 
them the so-called OPP technique (I18M21I) that allows for a fully numerical 
evaluation of the one-loop virtual amplitude. 

In this paper we describe one of the computational frameworks emerged 
during the last years, namely the HELAC-NL0. It incorporates several pieces of 
developed software, including HELAC-PHEGAS (043), CutTools (0), HELAC- 
DIP0LES (25), OneLOop (26), that have been already public for some time and 
HELAC-1L00P that is presented in this paper. This program in its current 
form is capable to evaluate fully numerically virtual QCD corrections to 
scattering amplitudes involving up to 7 particles directly attached to the loop 
composed by strongly interacting particles (gluons and quarks). In section [2] 
we will briefly describe the underlying algorithm providing information on the 
structure of the code that may be useful for potential developers. In the next 
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section [21 detailed instructions on running the code are given. Benchmark 
results can be found in section HJ along with a description of HELAC-DIPOLES 
in order the reader and potential user to have a complete overview of the 
full software. Finally in section [5] we critically review the current level of 
sophistication and the potential improvements. 



2. The HELAC-1L00P algorithm 

The aim of the program is to numerically evaluate the virtual contribu- 
tions needed in a next-to-leading order calculation. To this end the evaluation 
of both tree-order and one-loop amplitudes is necessary. 



The tree-order calculation algorithm is described in details in (122H24J) . 
In HELAC-1L00P the first step is the construction of the so-called tree-order 
skeleton, that contains all information for the evaluation of the amplitude. 
The primary input to this construction is the flavor (ifl(l:n)) assignment 
of the n (n) external particles. A special file manageable by the user, also 
used in HELAC, named constants. h, is used to numerically fix all physical 
constants needed (physics.f). 

The first action (helac_init in mastef _new. f ) is to enumerate (ncc) and 
define all possible color connections (icol (1 :n, 1 : 2) ). For a typical process 
consisting of n q numbers of (outgoing) incoming (anti) quarks and n g gluons, 
the number of color connections is a priori set to (n q +n g )\. Based on these 
data, the program is now constructing (within panl.f) using a top-down 
approach, all currents^ (list ( 1 : ngues , 1 : 18 , 1 : ncc) ) needed in a Dyson- 
Schwinger recursive representation of the amplitude, using the appropriate 
vertex functions (v3 for 3-vector-boson vertex, v4 for 4-vector-boson vertex, 
vff for vector-boson-fermion vertex, etc.), as shown schematically in the 
following figure 



4 Notice that the number of the needed currents is not a priori known; ngues. h is 
containing a user-defined estimate of it. 
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+ 




In In 

Notice that if, in partitioning the particles within blobs in the above 
diagram we keep the order of particles untouched, this is nothing but the 
so-called Berends-Giele (1271 ) recursive representation for color ordered am- 
plitudes. At the end of this skeleton construction the set of all currents 
for all color connections is stored in list. Notice that the number of color 
connections (ncc) can now be less than its a priori defined value. 

The one-loop n— particle amplitude can schematically be decomposed in 
a sum over terms of the form (m = 1, . . . , n) 



^~ d d d q N(q) 
with (i-dimensional denominators 



(1) 



A(?) = (q + Pi 



(2) 



where q is the loop momentum in d dimensions and N(q) is the numerator 
calculated also in d dimensions 0. The sum includes of course all terms with 
different loop-assignment structure: two structures may differ either by the 
number of denominators or by the different flavor and momenta appearing 
in the denominators. In that sense a closed gluon, ghost or massless quark 
loop, for instance, with the same momentum flow, is considered as different 
structure, although the denominators are identical. For the highest number of 
denominators each loop-assignment structure (taken into account the flavor 



5 When speaking about numerator function, it should be kept in mind that it generally 
contains propagator denominators not depending on the loop momentum 
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of the particles running in the loop) corresponds to a unique Feynman graph , 
but for m < n a collection of Feynman graphs with common loop-assignment 
structure should be understood. 

It is a well known fact that when d — >■ 4 limit is taken, the amplitude can 
be cast into the form 

A = E di Boxj + d Triangle^ + / J b i Bubble^ + E a * Tadpolej + R , (3) 

i i i i 

where Box, Triangle, Bubble and Tadpole refer to the well known scalar 
one-loop functions and R = R\ + R2 is the so-called rational term. 

The reduction of Eq.fJTJ to Eq. (JSJ) is the first ingredient of any approach 
aiming at the calculation of virtual corrections. In the following we will 
follow the so called reduction at the integrand level, developed by Ossola, 



Papadopoulos and Pittau (1181). The main idea is that any numerator function 
can be written as 



m—1 m—1 

N (q) = E d(i iii 2 i?) + d(q; ioi^iz) Yl A 

m—1 m—1 

+ [c(io*i*2) + c{q; i iii 2 )] JJ A 

m—1 m—1 

+ E [&(*o*i)+6(9;io*i)] II A 

m—1 m—1 

+ J2 t a (*o) + a(q;i )} \[ A 

m—1 

+ P(q) II A • (4) 

where now iV(g) and Di(q) are the four-dimensional versions of N(q) and 
Di(q). The coefficients d, c, 6 and a appearing in Eq.flU) are independent of 
the loop momentum and the same as the ones in Eq.([3]), whereas the new 
coefficients d, c, b, a and P(q), called also spurious terms, are depending on 
the loop momentum and they integrate to zero. 

Depending on the reduction method used, the calculation of any one-loop 
amplitude is placed in a very different perspective. For instance Eq.fllJ) can 
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be solved by computing the numerator functions for specific values of the 
loop momentum, that are solutions of equations of the form 



Di(q) = 0, for % = 0, . . . , M - 1 (5) 

It is customary to refer to these equations as quadruple (M = 4), triple 
(M = 3), double (M = 2) and single (M = 1) cuts. 

Calculating the numerator function for specific values of the loop momen- 
tum, opens the possibility to use tree-level amplitudes as building blocks. The 
reason is rather obvious: the numerator function is nothing but a sum of in- 
dividual Feynman graphs with the given loop-assignment structure and as 
we will see in a while, it is part of a tree amplitude with n + 2 particles. 
This is by itself a very attractive possibility, since one can use existing algo- 
rithms and tools that perform tree-order amplitude calculations, exploiting 
their automation, simplicity and speed. Indeed in the sequel we will describe 
how using HELAC we can also compute any one-loop amplitude. 



For the one-loop amplitude a skeleton construction is also performed (1281 ) . 
For a given external configuration (n, if 1, icol) the construction of all 
topologically inequivalent partitions (i.e. permutations) of the external par- 
ticles into all possible number of sets (blobs), is performed (loop.f and 
loop/loop_new. f 90). One such contribution is schematically represented in 
the following figure 




For those familiar with HELAC the numbering of external particles follows 
the binary representation used. In the present example 6 particles in direct 
contact with the loop are considered (hexagonal topologies). The allowed 
particle flavors (flavors. h) and colors running in the loop are defined (sub- 
routines check7, check6, etc.): the labels f and c in the figure above refer 
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to the possible flavor and color of the internal loop-particles. This construc- 
tion will continue to include also pentagon topologies, tetragonal topologies, 
triangle topologies, and bubble topologies. A typical collection of possible 
contributions, looks like 




Concerning the loop-momentum flow in these constructions, the conven- 
tion we have chosen is that it runs counterclockwise, and the loop-propagator 
connecting the blob that includes the particle number 1 and the last blob, is 
identified as D (q) of Eq.flT]). 

The selection of all the above mentioned contributions is enough for the 
calculation of the one-loop amplitude. To help the reader to understand the 
concept, the construction we have followed is equivalent to draw all possible 
one-loop Feynman graphs, and then collect them in sub-classes that are char- 
acterized by a common loop-assignment structure (after possible momentum 
shifts). 

In practice now, each numerator contribution, will be calculated as part 
of the n + 2 tree-order amplitude subject to the constraint that the attached 
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blobs will contain no propagator depending on the loop momentum 





and no denominator will be used for the internal loop propagators. Cutting 
now the line connecting the blob containing the particle number 1 and the 
last blob, it is easy to see that we have nothing more that a part of the 
n + 2 amplitude. The 'cut' particles, with flavor f and color connection 
(icol) appropriately defined (1281 ) . will now acquire their usual numbering of 
external particles in HELAC, namely 2 n and 2 n+1 , (64 and 128 for n = 6). 





8 
16 
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16 
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Using the above data the program knows how to reconstruct (subroutines 
numX_cuttools, X= 2 ... 7, in loop.f ) all information needed for the calcu- 
lation and store it as a sequence of sub-amplitudes or currents (listnum), 
quite as in the tree-order calculation, the main difference being that for 
a given color connection we have generically more than one contributions 
(nonums) characterized by different partitions of external particles, as well as 
flavors and colors running in the loop. 

Rational terms are classified in two categories: R\ are evaluated by 
CutTools (cts_xcut : rati); the R2 rational terms are calculated through 
Feynman rules as established in (|29h32|). In HELAC- 1L00P: physics. f in 
addition to all Standard Model couplings defined as in HELAC, couplings re- 
lated to the i?2 rational terms are also defined. Although we have to deal 
with a tree-order construction, the skeleton is build up following a procedure 
similar to that of the loop amplitude. The reason is that the special R2 
vertex has to appear only once for arbitrary scattering process at one loop. 
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Moreover these R2 vertices can join up to 4 particles. Therefore the skeleton 
construction (loop.f) starts with the distribution of all external particles in 
4, 3 and 2 sub-sets (blobs). Then these blobs, that represent currents or sub- 
amplitudes, are recursively (loop.f :mumX_r2 with X=2,3,4) defined in terms 
of the external particles and SM vertices, using the basic HELAC algorithm. 
The numerical evaluation of these contributions is identical to that of the 
tree-order amplitude. 

The data collected so far are used to evaluate numerically the amplitude 
in the so-called second stage. At this moment 4-momenta have to be sup- 
plied. In HELAC- 1L00P we provide three ways for dealing with 4-momenta: 
the testing modes (irambo=l, iint=l) in which case the 4-momenta are gen- 
erated randomly and one can use them for a number (nmc) of evaluation of 
matrix element and (irambo=0) in which case the 4-momenta are provided 
by the user through a file (mom); the third way (irambo=2, iint=0), the so- 
called re-weighting mode, is working through the reading of a Les Houches 
Events (LHE) file generated by the HELAC itself. 

Once 4-momenta are supplied, the external wave function vectors (4- 
dimensional complex vectors from wave.f) are computed. Then the tree- 
order amplitude for each color configuration is evaluated (panl .f :nextq). 
The total matrix element squared is then evaluated (master_new.f ) accord- 
ing to the standard formula 

J2KC a ^A a ,. (6) 

The color matrix C (rmatrix) has been already computed and stored in the 
first stage. 

For the one-loop amplitudes we follow the OPP reduction (nextql). To 
this end the CutTools program is interfaced and used (cuttools . intf . f ). 
The input to CutTools is the numerator function (numerator) evaluated 
at values of the loop momentum provided by CutTools itself. In order to 
compute the numerator function, the polarization vectors of the two extra 
'external' particles, after the one-particle cut, are also calculated. Within the 
Feynman gauge for gauge bosons, the sum over four different (4-dimensional) 
polarizations, that satisfy YU e i et i = 9^ ls performed. Ghost particles are 
also included. Finally for fermions, four vectors in spinor space, satisfying 
YiuQup — {i + m) a i3 are used. 

In the re-weighting mode, the actual calculation of the virtual corrections 



is organized using a re- weighting technique (1331 . |34J). To explain how this 
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works, let us start with the following equation (M. is the tree-order and £, 
the one-loop virtual matrix elements) 

°ab° +V = J dxrfxzd^UxJMxz) (\M\ 2 + MC* + M*C) , (7) 

which gives the sum of leading order (LO) and virtual (V) contributions for 
a scattering ab — > m-particles. It can be re-written as 

r f SAC* 4- SA*C\ 
°a b ° +V = J dx 1 dx 2 d$ m f a (x 1 )f b (x 2 )\M\ 2 f 1 + J . (8) 

Since £ is a time consuming function one would like to calculate it as few 
times as possible. To this end a sample of un-weighted events is produced 
based on the tree order distribution, namely 

1 da LO 

9{X) = 9ixi , X2 ,* m) = — 1 -^-, (9) 

satisfying / dXg(X) = 1. The sample S of un-weighted events has the fol- 
lowing property, 

-L^O^H fdXg(X)0(X), (10) 
i€S J 

where the equality should be understood in the statistical sense, and 0(X) 
is any well-defined function over the integration space. Now it is trivial to 
see that if 



then 



1 a LO+v 



atXWH-As-- (12) 



N S i&S ' a ab 



In practice the sample of tree order un-weighted events includes all infor- 
mation on the integration space, namely, the color assignment, the (random) 
helicity configuration, the fractions X\ and x 2 an d the m— body phase-space. 
For future convenience it is produced in a standard Les Houches format (j3oT ). 
One-loop contributions are only calculated for this sample of un-weighted 
events, and the weight assigned to each of those events is given by 

MC* + M*C 
W= \M\ 2 ■ (13) 
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The total virtual contribution can now be easily estimated by 



a v = (w) a B , (14) 

where a B is the Born cross section, already included in the LHE file. More- 
over, the sample of events including the information on w, can be used to 
produce any kinematical distribution, according to the Eq. [TOJ 

Most reduction algorithms suffer from numerical instabilities usually cau- 
sed by the presence of small Gram determinants. In order to detect the phase- 
space points with unstable behavior with respect to the reduction procedure, 



one can rely on several tests (36). In the current version of HELAC-1L00P 
we employ actually the so-called gauge test. Our experience shows that 
the effect of numerical instabilities is more pronounced when higher rank 
tensor integrals are involved (namely high powers of loop momentum in the 
numerator function). This is often correlated with the presence of gluons in 
the scattering amplitude under consideration. In that case the gauge test 
can be performed and it is equivalent to replace the polarization vector of 
an external gluon with its momentum. Events that fail to obey the Ward 
identity are separated from the initial sample. They can be treated then in 
higher (quadruple) precision. In any case experience shows that only for very 
complicated processes, like for instance gg — > tigg, the effect is appreciable. 
Care should be taken in properly defining the numerical criterion of rejection. 
For the moment we replace the gluon polarization vector with the gluon 
momentum normalized to its energy e M — > p^/p and reject events when the 
computed matrix element differs from zero by an amount greater to 10~ 9 . 
On the other hand we should emphasize that this procedure cannot be the 
final answer, since it has a limited applicability depending on the presence 
of external gluons in the scattering process under consideration. 

We should also notice that for being able to test the calculations we 
provide the numerical evaluation of the infraredpart (wl_I) of the so-called 



/-operator. This is achieved as described in ( 12 8l ) by computing the color- 



correlated matrix (rmatrix_l) and then using the tree-order amplitudes. 



3. How to run the code 

The code is written in Fortran 90 and needs no additional software to 
run. We use the fortran compiler gf ortran as the default, but in many ap- 
plications we have used also lahey95 and if ort. Unpacking the distributed 
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tarball helaclloop.tgz will create the directory HELAC1L_0FFICIAL. In this 
you will find the subdirectories examples, run, src/lLOOP, src/utils and 
src/TREE. For the convenience of the potential user, three script files named 
run_testing, run_reweight and run_GC can be found in the subdirectory 
run that can guide the user to run the code as described in details below. 

3.1. 1LOOP 

• Edit the file constants. h and define your own physical parameters if 



needed (see (I24l ) for explanation). 



Run the script file scriptmake 

./scriptmake nl n2 n3 n4 where nl means normal double preci- 
sion and 1 quadruple (if supported by fortran version), n2 for compil- 
ing everything from scratch and 1 only the latest files, n3 is empty by 
default and can be used to pass additional flags to fortran commands 
(see the scriptmake file), and finally n4 either empty or GC that means 
gauge check mode. For instance the command 



./scriptmake "" "GC" 

will produce the main_onep_dpGC . exe executable file, for use in gauge 
check, where 

./scriptmake 

will produce the main_onep_dp . exe executable file. 

Edit the input file input and define appropriately the input parameters: 

— iint for re-weighting mode, 1 for all others 

— ibv for full summation over colors, 1 for Monte Carlo over colors 

— iverbose 0...3 different levels of verbosity 

— repeat 1 only skeleton construction, 2 only amplitude evaluation 
(assuming the skeleton is present), for one shot 

— iranhel for full summation over helicities 1 for MC over helici- 
ties 

— n number of particles 
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— flavors flavor of particles according to the HELAC list 

— if lag for internal definition of helicities 1 for user providing 
helicities 

— ihiggs if no Higgs is included 1 if Higgs is included if allowed 

— loopi T if loop amplitude is calculated F if only tree is calculated 

— only T if only QCD couplings are allowed, F if also EW couplings 
are allowed 

— withqcd T if QCD is included, F elsewhere 

— irambo if momenta are provided by the user (through mom file), 
1 if are generated randomly by RAMBO, 2 if provided by the .lhe 
LHE file generated by HELAC 

— e energy in GeV 

— nmc number of phase-space points to be evaluated (for irambo 
or 1) 

— mom file with the 4-momenta in the format E, p x , p y , p z , m 

— momout output file for user provided momenta 

— muse ale the renormalization scale 

3.2. TREE 

In this directory you can generate the LHE file, for later use in the re- 
weighting procedure. The generation is identical to the standard HELAC-PHEGAS 



procedure. So we refer the potential user to (24). Some new elements with 
respect to the standard treatment have been added, without affecting the 
generation procedure. 

• In this version a color MC is being used according to (5). There is a 
new keyword color_f lag for a full color summation and 1 for a MC 
over colors. 



An interface to phase-space generator KALEU (37) is also present, the 



user can set the phasespace_f lag to for PHEGAS or 1 for KALEU. 

• Finally the oneloop_rewgt keyword has to be set to true (T) in order 
the generated LHE to be used for the re-weighting procedure described 
above. 

The tree-order generation will result to a LHE file, sampleGOGOTqTa. lhe 
for instance for the process gg — > ti. 
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4. Results and Benchmarks 



In this section we will give characteristic examples of running the code. 
The tarball file also contains the examples directory where these results are 
stored. 

4-1. The testing mode 

Let us choose a relatively simple example ud — > W + + ng with n — 1, 2, 3. 
We start with simplest case n — 1. The input file looks like: The output 
file (the reader can find it as uDWg_gen in the official distribution examples) 
incorporates information on how the different currents are formatted. The 
main result of this run is encoded in the files 

• tree_UqDaW+GO . in including all the information needed by HELAC-1L00P 
to evaluate numerically the amplitude and 

• treeli_UqUaW+GO . in with some brief description on the virtual ampli- 
tude generation 

Now by changing the repeat keyword value from 1 to 2, we have the 
numerical evaluation of the amplitude in uDWg_out. Most of the content of 
the output file is self-explanatory. We focus here on its major aspects. After 
the printout of physical constant the user (if iverbose > 1) will see the 
following 

UqDaW+GO 

a line printing out the process under consideration ud — > W + g. Then a 
bunch of lines like 



INFO ============================================= 

INFO COLOR 1 out of 2 

INFO 26 -3 51 14 34 32 -4 2000 1 10 

INFO 2 10 -4 61 18 35 42 -4 2000 1 13 

INFO 2 10 -4 6018 35 42 -4 2000213 

INFO 2 14 -3 7 1 2 4 34 3 10 -4 6 1 1 

INFO 2 14 -3 7228 35 46 -3 50001 13 

INFO 2 14 -3 7028 35 46 -3 5000213 

where each line represents a sub-amplitude needed in the construction of 
the tree-order amplitudes corresponding to the 1 out of 2 color 

connection configurations existing for this process. Then after exhausting all 
tree-order information you start seeing 



15 



LOOP T 

INFO ==================================== 

INFO COLOR 1 out of 2 



INFO 


number 


of nums 






33 






















INFO 


NUM 






1 


of 






33 






6 














INFO 


3 


24 


3 


7 


1 


1 


8 


35 


4 


16 


3 


5 











1 


1 


INFO 


3 


24 


3 


7 





1 


8 


35 


4 


16 


3 


5 











2 


1 


INFO 


3 


28 


4 


8 


1 


1 


4 


34 


3 


24 


3 


7 











1 


1 


INFO 


1 


30 


35 


9 


1 


1 


2 


-4 


2 


28 


4 


8 














1 


INFO 


2 


62 


-3 


10 


1 


1 


30 


35 


9 


32 


-3 


6 











1 


1 


INFO 


2 


62 


-3 


10 





1 


30 


35 


9 


32 


-3 


6 











2 


1 


INFO 


4 


8 


4 


2 


1 


3 


3 


4 


35 























3 



INFOYY 1 

again for the 1st color connection. As it is printed there are now 33 
(number of nums 33) contributions and the first one is nothing 

but the following box graph 




The relevant topological information can easily be read off from the line 
INFO 48421334 35 000000031 

the first number after the keyword INFO, namely 4, being the number 
of loop propagators, then 8 4 2 1 the ordering of blobs and finally 
3 3 4 35 (uudg) the corresponding flavor assignment. 
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Out of the 33 contributions, the last 6 refer to R 2 tree-order like ones 
including a unique special vertex each. In the following example a special R 2 
udW vertex 



INFO NUM 29 of 33 3 
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INFOYY 1 

is used as shown also schematically in the following Feynman graph: 




At the end of the file, after momenta have been generated and printed, 
you find the most important result which is the numerical values of matrix 
elements: 

HELICITY CFGS 6 

total amplitude squared LO = 0.73261294162118751 

total amplitude squared U0 = 4 . 88872424257442914E-002 

total amplitude squared TO = 4 . 88872424257442914E-002 
ratio = 6.67299738352457872E-002 

total amplitude squared Ul = 2 . 38831983483179001E-002 

total amplitude squared Tl = -2 . 26466317876591848E-002 

total amplitude squared T2 = -6 . 87832271575313514E-002 

total amplitude squared II = -2 . 26466317876590842E-002 
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total amplitude squared 12 = -6 . 87832271575313653E-002 
total amplitude squared RO = 0.0000000000000000 
total amplitude squared Rl = -4 . 65298301359770849E-002 

In this case there are 6 non-zero helicity configurations, L0 is the lead- 
ing order \Ai\ 2 properly summed and averaged, U0 is the finite part of the 
one-loop amplitude not including coupling constant renormalization (but in- 
cluding mass renormalization if present), TO the renormalized one, and so 
on for Ul, Tl for the - and U2, T2 for the 4 terms. The II and 12 are 
the corresponding poles predicted by the /-operator and the level of agree- 
ment between T2 and 12 as well as between Tl and II reflect the precision 
achieved. 

In the website of HELAC-NL0, results can also be found for ud — > W + gg 
and ud — > W + ggg. 

4-2. Re-weighting mode 

The distributed version contains reproducible results for the process gg — > 
ti. To obtain those results you have to first generate in TREE sub-directory 
the LHE file by just executing the following command 

./run.sh user.inp myenv 

as in the usual HELAC-PHEGAS. In general the user has to edit the user . inp 
input file and define the corresponding input parameters. 

Copy the produced sample file sampleGOGOTqTa. lhe in the same di- 
rectory where the main_onep_dpGC . exe is used in case you have to run a 
gauge test. To this end the set-up of the input parameters has to be as 
follows: iint= 0, ibv= 1 and irambo= 2. When now the executable is run- 
ning two sample LHE files will be generated sampleGOGOTqTa.GAUGECHECK- 
_FAILED . lhe and sampleG0G0TqTa_GAUGECHECK_PASSED . lhe, containing the 
events that have failed and passed the test respectively. 

Then rename this last file to sampleGOGOTqTa. lhe and run 

. /main_onep_dp . exe < input ><<your output» 

The final result is the LHE file sampleG0G0TqTa_WEIGHTED. lhe which can 
now be used in conjunction with the HELAC-DIP0LES to get the full NLO 
corrections. 
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The events if any, in sampleGOGOTqTa_GAUGECHECK_FAILED . lhe 

can be reprocessed in higher numerical precision^]: 

. /main_onep_qp . exe < input ><<your output» 

The resulting file is again sampleGOGOTqTa_WEIGHTED . lhe which can now 
be combined!! with the one produced with the passed events. Care of course 
has to be paid by the user in properly managing and re-naming the emerging 
files. 

4.3. EEL AC-DIP OLES 

In this subsection we will briefly review the use of HELAC-DIPOLES in order 
the user to be able to have a more complete overview of the full software. 



For a detailed description please refer to (1251 ) . 

Unpack the tarball dipoles.tgz inside the main HELAC directory and 
copy alphas_std.h from the dipoles directory to HELAC. This sets up two- 
loop running consistently with the PDFs used by HELAC-DIPOLES. 

• Configure HELAC -PHEGAS 

— Edit file myenv to setup the libraries and compilers, 

— Edit user . inp, to specify the input physical parameters, 

— Edit getqcdscale .h, if a running coupling constant is needed, 

— Edit nh . h, to specify the number of histograms that will be gen- 
erated. 

Apart from myenv, for all other files described so far the user may use 
their default configuration. To run 

./run.sh user. inp myenv 

This procedure will set-up several parameters needed by HELAC-DIPOLES 
later on, among others the cuts to be used by the integrated dipoles 
(to be run with make run_I and make run_KP, see below). 



Such a possibility can be realized with lahey95 and ifort that incorporate at the 
compilation level the option of quadruple precision. See also the script file scriptmake 
for more details. 

7 In the sub-directory src/utils/COMBINE useful tools for combining LHE files are avail- 
able. 
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• Configure HELAC-DIPOLES 



— Edit dipoles . input, to specify the process and various optimiza- 
tion parameters just as in the original HELAC, 

— Edit dipoles . conf, to set specific parameters for the calculation 
of the dipoles, as follows: 

* onlyreal: if set to true, only real corrections for a given pro- 
cess will be calculated. The cuts are then specified in cuts.h. 
The result must coincide with the one of the original HELAC 
with the same input parameters. This option is included for 
testing purposes. 

* onlylast: if set to true, only those dipoles will be included, 
which contain the last particle (for correctness it must be a 
parton). This is useful for some processes, where it is clear 
that only the last particle can be soft / collinear, and the book- 
keeping remains simple to obtain the full result at NLO. 

* onlydiv: if set to true, only divergent dipoles will be included. 
Non-divergent dipoles correspond to a pair of massive quarks 
in the final state. They are only useful to get rid of large 
Sudakov logarithms, but are not essential for the finiteness of 
the real radiation contribution. 

* hybrid: if set to true, non-parton polarizations will be summed 
over by a continuous Monte Carlo integration over a phase pa- 
rameter. 

* signmode: defines how positive and negative contributions 
are to be treated: - the result is left unchanged, whether 
positive or negative, 1 - only positive numbers, a negative 
result is set to zero, 2 - only positive numbers, but the sign 
of a negative result is changed, and positive results are set to 
zero 

* sumtype: in the first phase (preferably for phase space opti- 
mization), the summation over helicities of the partons can 
be performed in three different ways: - exact fast summa- 
tion (independently for real radiation and dipoles), 1 - ex- 
act slow summation (for a given helicity configuration both 
real radiation and the dipole sum will be calculated), 2 - flat 



20 



Monte Carlo summation over all non-vanishing helicity con- 
figurations. In practice, it is recommended to use the last 
option. 

* nsumpol: number of accepted points to be summed over he- 
licity with the method specified by sumtype. The counting 
starts after phase space optimization is finished. 

* noptpol: number of accepted points to be used for helicity 
sampling optimization. During helicity sampling optimiza- 
tion, slow summation over helicity configurations (in the sense 
defined in the description of sumtype) is performed. It is 
therefore recommended to keep this number relatively small 
(of the order of a few hundred to a thousand). 

* nuptpol: number of accepted points after which an updated 
of the helicity sampling weights is performed. This number 
should be rather large for best results (at least an order of 
magnitude larger than noptpol). 

* alphaMinCut: lowest value of alphaMin, below which a point 
will be rejected altogether, because of risk of numerical insta- 
bilities. For exact definition of alpha, see (jjjl). 




* alphaMaxII, alphaMaxIF, alphaMaxFI, alphaMaxFF, kappa: 



— j etf unctions . f , for non-trivial jet functions, although most work 
should be performed on FJmpo (jet function for real radiation), 

— cuts.h, to specify general cuts to be used by the jet functions, 

— histograms . f, to define histograms, 

— seed, input, to change the random number seed, which is useful 
for trivial parallelization runs (this should be an integer). 

Again one can use the existing default configuration. To run the pro- 
gram 

export FC="<<your fortran>>" 
make 

make run 

This will compile the dipole subtracted version and run it. In case this 
is the first run for a given process, the calculation will be stopped and 
you should run 
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make trees 
make run 



The former will generate and store in the subdirectory helac_trees 
all skeleton files for the different subprocesses needed by the dipole- 
subtraction. 

For the / operator the corresponding input and configuration files are 
dipoles_I . input and dipoles_I . conf . The parameters defined in 
the configuration file are: the renormalization scale, the number of light 



Nf and heavy Np quark flavors, used in the definition of / operator (125 



381 ) . the a parameter that controls the integration of dipole functions 
over the available phase-space and the parameter k defined in (j38l ). To 



run 



make run_I 



• For the K + P operator the corresponding input and configuration files 
are dipoles_KP . input and dipoles_KP . conf. In addition to the com- 
mon parameters described so far for the / operator, the configuration 
file contains also the definition of the factorization scale and the num- 
ber and flavors of initial state partons to be taken into account. To 
run 



make run_KP 



The result of the run consists of several files. A typical output will have 
at the end the following information 

out of 1000000 1000001 points have been used 

and 340079 points resulted to =/= weight 

whereas 659922 points to weight 

estimator x: 0.381485D-16 

estimator y: 0.143375D-32 

estimator z: 0.205540D-65 

average estimate : 0.381485D-16 
+\- 0.378649D-16 

variance estimate: 0.143375D-32 
+\- 0.143367D-32 

be aware that the error estimate may be bad! 

estimator x: 0.835522D-04 
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estimator y: 0. 156118D-12 
estimator z: 0.930801D-29 
average estimate : 0.835522D-04 
+\- 0.395118D-06 
variance estimate: 0.156118D-12 
+\- 0.305090D-14 
total XS -8.35522401509068940E-005 3 . 95117621007448708E-007 
lwri: points have used 0.0000000000000000 

2212 2212 7000.0000000000000 7000.0000000000000 3 

I error: 99.256586045324298 
7„ error: 0.47289889570121785 

which states that from a run of 1 million points, 340,079 have been used 
after cuts and the positive part of the cross section is 0.381485D-16 with 
a statistical uncertainty +\- 0.378649D-16 whereas the negative part 
is 0.835522D-04, +\- 0.395118D-06 (in nanobarns). See also Q for 
more explanations. Also files named hi_f ile, hi_f ile_I, hi_f ile_KP will 
be generated with all data needed for the histograms. 



5. Outlook 



The progress in NLO calculations seen over the last years has made 
the development of an automatized computational framework a realistic 



leveiopme 

task (0,H0,|2^,l39H42|). In as much as in the tree-order generation NLO 
programs will be able to generate LHE files ready for use in physics analyses. 
Nevertheless there are several open issues to be addressed in the near future: 

• For the moment LHE can be generated for tree-order plus virtual cor- 
rections. For real corrections usually one has to deal with so-called 
"weighted" events, namely a large collection of phase-space integration 
points. Taking into account also that the Monte-Carlo convergence of 
the real corrections, especially the real-subtracted part is quite slow, a 
solution will be very welcome. Within HELAC-NL0 we plan to further 
investigate the possibility of using alternative phase-space algorithms, 
sampling over colors and other subtraction methods in order to achieve 
a significant improvement in the overall efficiency. 



The interface to Parton-Shower (143M46I) programs is also an impor- 
tant issue at the phenomenological level. There are already several 
steps taken towards this direction. So far, HELAC-1L00P has been in- 
terfaced to the P0WHEG-B0X framework at the purpose of studying spe- 
cific processes like pp — > ttH ( 1471 ) and pp — > ttj (1481). including NLO 
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QCD corrections matched to a Parton Shower evolution, followed by 
hadronization and hadron decay, up to final predictions at the hadron 
level to be compared to LHC and Tevatron data. We plan to further 
investigate this issue with the aim to integrate and automatize the full 
procedure. 

• To address the incorporation of the full set of Electroweak corrections 
as well as theories beyond the SM, a re-implementation of the HELAC 
algorithm is desirable. In such a process, we aim also to include several 
straightforward improvements in constraining the redundancy of the 
actual computation, resulting to a significant reduction in computing 
time and resources. 

Finally it will be very interesting to advance beyond one loop. Both the 
OPP reduction method and the recursive approach to scattering amplitude 
computation may open the road to highly efficient calculations at the two- 



loop level (1491). 



Acknowledgments: It is our pleasure to acknowledge very fruitful dis- 
cussions with Z. Bern, D. Kosower, L. Dixon, H. Ita, D. Maitre, Z. Kun- 
szt, A. Denner, C. Anastasiou, G. Ossola, P. Mastrolia, W. Giele, R. Fred- 
erix, A. Lazopoulos, P. Draggiotis and I. Malamos. We also acknowledge 
partial support from FALCON FPA 2008-02984, HEPTOOLS MRTN-CT- 
2006-035505, LHCPhenoNet PITN-GA-2010-264564. M.C was supported by 
the Heisenberg and by the Gottfried Wilhelm Leibniz Programmes of the 
Deutsche Forschungsgemeinschaft. M.W acknowledges support by the Ini- 
tiative and Networking Fund of the Helmholtz Association, contract HA-101 
(Physics at the Terascale). M.V.G acknowledges the University of Granada, 
Departamento de Fisica Teorica y del Cosmos, and the NCSR Demokri- 
tos, Institute of Nuclear Physics, for hospitality. C.G.P and R.P would like 
to thank the KITP at UCSB for the kind hospitality offered during the 
completion of this paper. A. K. was supported by the TAMOP 4.2.1./B- 
09/1/KONV-2010-0007 and the TAMOP 4.2.2/B-10/1-2010-0024 projects. 

References 



24 



References 



[1] J. R. Andersen, et al., The SM and NLO multileg working group: Sum- 
mary reporl larXiv : 1003 . 12411 

[2] C. F. Berger, et al., An Automated Implementation of On-Shell 
Methods for One- Loop Amplitudes, Phys. Rev. D78 (2008) 036003. 



arXiv : 0803 .4180, doi : 10 . 1103/PhysRevD . 78 . 036003 . 



[3] R. K. Ellis, W. T. Giele, Z. Kunszt, K. Melnikov, G. Zanderighi, One- 
loop amplitudes for W+3 jet production in hadron collisions, JHEP 01 



(2009) 12. [arXiviOBl 0.2762, doi : 10 . 1088/1126-67 08/2009/01 /012 



[4] C. F. Berger, et al., Precise Predictions for W+3 Jet Produc- 
tion at Hadron Colliders, Phys. Rev. Lett. 102 (2009) 222001. 



larXiv: 0902 .27601 |doi : 10 . 110"37 PhysRe vLett . 102 . 222001 [ 



[5] G. Bevilacqua, M. Czakon, C. G. Papadopoulos, R. Pittau, M. Worek, 
Assault on the NLO Wishlist: pp -> tt bb, JHEP 09 (2009) 109. 
larXiv : 0907 . 4723| [doi : 10 . 1088/1126-6708/2009/09/109] 

[6] G. Bevilacqua, M. Czakon, C. G. Papadopoulos, M. Worek, Dom- 
inant QCD Backgrounds in Higgs Boson Analyses at the LHC: 
A Study of pp -> t anti-t + 2 jets at Next- To-Leading Or- 
der, Phys. Rev. Lett. 104 (2010) 162002. larXiv: 1002.40091 
|doi : 10 . 1103/PhysRevLett . 104. 162002 [ 

[7] K. Melnikov, M. Schulze, NLO QCD corrections to top quark 
pair production in association with one hard jet at hadron col- 
liders, Nucl. Phys. B840 (2010) 129-159. larXiv: 1004~3 284. 
|doi : 10 . 1016/ j . nuclphysb . 2010 . 07 . 003[ 

[8] A. Bredenstein, A. Denner, S. Dittmaier, S. Pozzorini, NLO 
QCD corrections to pp — > t anti-t b anti-b + X at the 
LHC, Phys. Rev. Lett. 103 (2009) 012002. larXiv: 0905 .01101 
doi : 10 . 1103/PhysRevLett . 103 . 012002| 



[9] A. Bredenstein, A. Denner, S. Dittmaier, S. Pozzorini, NLO QCD cor- 
rections to top anti-top bottom anti-bottom production at the LHC: 
2. full hadronic results, JHEP 03 (2010) 021. larXiv: 1001 .40061 



doi : 10 . 1007/ JHEP03 (2010) 021 , 



25 



[10] G. Bevilacqua, M. Czakon, A. van Hameren, C. G. Papadopoulos, 
M. Worek, Complete off-shell effects in top quark pair hadroproduc- 
tion with leptonic decay at next-to-leading order, JHEP 02 (2011) 083. 
larXiv : 1012 . 42301 Idoi : 10 ■ 1007/JHEP02(2011)083| 

[11] A. Denner, S. Dittmaier, S. Kallweit, S. Pozzorini, NLO 
QCD corrections to WWbb production at hadron collid- 
ers, Phys. Rev. Lett. 106 (2011) 052001. larXiv: 1012.39751 
|doi : 10 . 1103/PhysRevLett . 106 . 052001 [ 

[12] G. Bevilacqua, M. Czakon, C. G. Papadopoulos, M. Worek, Hadronic 
top-quark pair production in association with two jets at Next-to- 
Leading Order QCD arXiv: 1108.28511 

[13] C. F. Berger, et al., Precise Predictions for W + 4 Jet Production 
at the Large Hadron Collider, Phys. Rev. Lett. 106 (2011) 092001. 



arXiv: 1009 . 2338, doi : 10 . 1103/PhysRevLett . 106 . 092001 



[14] H. Ita, et al., Precise Predictions for Z + 4 Jets at Hadron Collid- 
er darXiv: 1108.22291 

[15] T. Melia, K. Melnikov, R. Rontsch, G. Zanderighi, NLO QCD cor- 
rections for W + W~ pair production in association with two jets at 
hadron colliders, Phys. Rev. D83 (2011) 114043. arXiv: 1104. 23271 



doi : 10 . 1103/PhysRevD . 83 . 114043 . 



[16] G. P. Salam, Perturbative QCD for the LHC, PoS ICHEP2010 (2010) 
556. larXiv: 1103. 13181 

[17] R. K. Ellis, Z. Kunszt, K. Melnikov, G. Zanderighi, One-loop calcu- 
lations in quantum field theory: from Feynman diagrams to unitarity 
cutE farXiv: 1105.43T91 

[18] G. Ossola, C. G. Papadopoulos, R. Pittau, Reducing full one- 
loop amplitudes to scalar integrals at the integrand level, 



Nucl. Phys. B763 (2007) 147-169. arXiv:hep-ph/0609007 



doi : 10 . 1016/ j . nuclphysb . 2006 . 11 . 012 



[19] G. Ossola, C. G. Papadopoulos, R. Pittau, Numerical Evaluation 
of Six-Photon Amplitudes, JHEP 07 (2007) 085. larXiv: 0704 .12711 



doi : 10 . 1088/11 26-6708/2007/07/085 . 



26 



[20] G. Ossola, C. G. Papadopoulos, R. Pittau, CutTools: a pro- 
gram implementing the OPP reduction method to compute one- 
loop amplitudes, JHEP 03 (2008) 042. larXiv: 0711 .35961 
doi : 10 . 1088/1126-6708/2008/03/042. 



[21] G. Ossola, C. G. Papadopoulos, R. Pittau, On the Rational Terms 
of the one-loop amplitudes, JHEP 05 (2008) 004. llarXiv: 0802. 18761 
|doi : 10 . 1088/1126-6708/2008/05/0041 

[22] A. Kanaki, C. G. Papadopoulos, HELAC: A package to compute elec- 
troweak helicity amplitudes, Comput. Phys. Commun. 132 (2000) 306- 
315. |arXiv : hep-ph/0002082} |doi : 10 . 1016/S0010-4655 (00) 0015lTx| 

[23] C. G. Papadopoulos, PHEGAS: A phase space generator for automatic 
cross- section computation, Comput. Phys. Commun. 137 (2001) 247- 
254. |arXiv : hep-ph/0007335| |doi : 10 . 1016/S0010-4655 (01) 00163^11 

[24] A. Cafarella, C. G. Papadopoulos, M. Worek, Helac-Phegas: a generator 
for all parton level processes, Comput. Phys. Commun. 180 (2009) 1941- 
1955. llarXiv: 0710 .24271 |doi : 10 . 1016/ j . cpc . 2009 . 04 . 023] 

[25] M. Czakon, C. G. Papadopoulos, M. Worek, Polariz- 
ing the Dipoles, JHEP 08 (2009) 085. larXiv :0905~0 883. 
|doi : 10 . 1088/1126-6708/2009/08/0851 

[26] A. van Hameren, OneLOop: for the evaluation of one-loop 
scalar functions, Comput. Phys. Commun. 182 (2011) 2427-2438. 
larXiv : 1007 . 47161 |doi : 10 ■ 1016/j . cpc . 2011 . 06 ■ 011 1 



[27] F. A. Berends, W. T. Giele, Recursive Calculations for 
Processes with n Gluons, Nucl. Phys. B306 (1988) 759. 
|doi: 10 .1016/0550-3213(88) 90442-71 

[28] A. van Hameren, C. G. Papadopoulos, R. Pittau, Automated 
one-loop calculations: a proof of concept, JHEP 09 (2009) 106. 
larXiv : 0903 . 46651 |doi : 10 ■ 1088/1126-6708/2009/09/106] 



[29] P. Draggiotis, M. V. Garzelli, C. G. Papadopoulos, R. Pittau, Feynman 
Rules for the Rational Part of the QCD 1-loop amplitudes, JHEP 04 



(2009) 072. [arXiy_i0903. 0356, doi : 10 . 1088/1126-6708/2009/04 /072 



27 



[30] M. V. Garzelli, I. Malamos, R. Pittau, Feynman rules for 
the rational part of the Electroweak 1- loop amplitudes, JHEP 
01 (2010) 040, [Erratum-ibid. 1010:097,2010]. larXiv: 0910 .31301 
Idoi : 10 . 1007/ JHEP01 (2010) 040] 

[31] M. V. Garzelli, I. Malamos, R. Pittau, Feynman rules for the ra- 
tional part of the Electroweak 1- loop amplitudes in the gauge 
and in the Unitary gauge, JHEP 01 (2011) 029. arXiv: 1009.43021 
doi : 10 . 1007/ JHEP01 (2011) 029] 



[32] M. V. Garzelli, I. Malamos, R2SM: a package for the analytic compu- 
tation of the R2 Rational terms in the Standard Model of the Elec- 
troweak interactions, Eur. Phys. J. C71 (2011) 1605. arXiv: 1010. 12481 



doi: 10 .1140/epjc/sl0052-011- 1605^6] 



[33] A. Lazopoulos, K. Melnikov, F. Petriello, QCD corrections to tri-boson 



production, Phys. Rev. D76 (2007) 014001. arXiv:hep-ph/0703273 



doi : 10 . 1103/PhysRevD . 76 . 01 4001 . 



[34] T. Binoth, G. Ossola, C. G. Papadopoulos, R. Pittau, NLO QCD correc- 
tions to tri-boson production, JHEP 06 (2008) 082. larXiv: 0804 .03501 
|doi : 10 . 1088/1126-6708/2008/06/0821 

[35] J. Alwall, et al., A standard format for Les Houches event files, Corn- 



put. Phys. Commun. 176 (2007) 300-304. arXiv:hep-ph/0609017 



doi: 10. 1016/ j . cpc . 2006 . 11 . 010 



[36] R. Pittau, Testing and improving the numerical accuracy of the 
NLO predictions, Comput. Phys. Commun. 181 (2010) 1941-1946. 



arXiv: 1006.3773, doi : 10 . 1016/ j ■ cpc . 2010 . 09 .013 



[37] A. van Hameren, Kaleu: a general-purpose parton-level phase space 
generatoi tarXiv : 1003 . 49531 

[38] S. Catani, S. Dittmaier, M. H. Seymour, Z. Trocsanyi, The dipole 
formalism for next-to-leading order QCD calculations with massive 
partons, Nucl. Phys. B627 (2002) 189-265. |arXiv : hep-ph/0201036] 
|doi:10.1016/S0550-3213(02)00098-6l 



28 



[39] T. Gleisberg, F. Krauss, Automating dipole subtraction for QCD NLO 
calculations, Eur. Phys. J. C53 (2008) 501-523. larXiv: 0709 .28811 
|doi:10.1140/epjc/sl0052-007-0495-0| 

[40] C. F. Berger, et al., Multi-jet cross sections at NLO with BlackHat and 
Sherp diTrXiv: 0905 .27551 

[41] R. Frederix, S. Frixione, F. Maltoni, T. Stelzer, Automation of next- 
to-leading order computations in QCD: the FKS subtraction, JHEP 10 
(2009) 003. larXiv : 0908 . 4272l[doi : 10 . 1088/1126-6708/2009/10/0031 

[42] V. Hirschi, et al., Automation of one- loop QCD corrections, JHEP 05 
(2011) 044. larXiv : 1103 . 062l[ Idol : 10 . 1007/JHEP05(2011)044| 

[43] S. Frixione, B. R. Webber, Matching NLO QCD computations and par- 
ton shower simulations, JHEP 06 (2002) 029. |arXiv : hep-"ph70 204244 . 

[44] P. Nason, A new method for combining NLO QCD with shower 



Monte Carlo algorithms, JHEP 11 (2004) 040. arXiv:hep-ph/0409146 , 
|doi : 10 . 1088/1126-6708/2004/11/0401 

[45] S. Frixione, P. Nason, C. Oleari, Matching NLO QCD computations with 
Parton Shower simulations: the POWHEG method, JHEP 11 (2007) 



070. larXiv :0709T2 092. do i : 10 . 1088/1126-6 708/20 07/11/ 070 . 



[46] S. Alioli, P. Nason, C. Oleari, E. Re, A general framework for 
implementing NLO calculations in shower Monte Carlo programs: 
the POWHEG BOX, JHEP 06 (2010) 043. larXiv: 1002 .~258il 



doi : 10 . 1007/ JHEP06 (2010) 043 



[47] M. V. Garzelli, A. Kardos, C. G. Papadopoulos, Z. Trocsanyi, Stan- 
dard Model Higgs boson production in association with a top anti-top 
pair at NLO with parton showering, Europhys. Lett. 96 (2011) 11001. 
larXiv: 1108 .03871 |doi : 10 . 1209/0295-5075/96/11001] 



[48] A. Kardos, C. Papadopoulos, Z. Trocsanyi, Top quark pair production 
in association with a jet with NLO parton showering ar Xiv: 1101 .2 672, 

[49] P. Mastrolia, G. Ossola, On the Integrand- Reduction Method for Two- 
Loop Scattering Amplitude darXiv : 1107. 6041 . 



29 



